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The pion-mass dependence of the nucleon mass witliin the covariant SU{2) baryon chiral per- 
turbation theory both without and with explicit A (1232) degrees of freedom up to order p^ is 
investigated. By fitting to lattice QCD data in 2 and 2+1 flavors from several collaborations, for 
pion masses M^ < 420 MeV, we obtain low energy constants of natural size and compatible with 
pion nucleon scattering data. Our results are consistent with the rather linear pion-mass dependence 
showed by lattice QCD. In the 2 flavor case we have also performed simultaneous fits to the nucleon 
mass and ct^jv data. As a result of our analysis, which encompasses the study of finite volume 
corrections and discretization effects, we report a value of ct^jv = 41(3) (1) MeV in the 2 flavor case 
and a-nN ~ 52(3) (8) MeV for 2+1 flavors, where the inclusion of the A (1232) resonance changes the 
results by around 9 MeV. In the 2 flavor case we are able to set independently the scale for IQCD 
fT^ . data, given by a Sommer scale of ro = 0.493(23) fm. 
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I. INTRODUCTION 



The nucleon mass A/jv is one of the fundamental observables in nature. It arises from the complex and not well 

r^ \ understood quark-gluon dynamics in the non-perturbative regime of quantum chromodynamics (QCD). Nevertheless, 

QiH' important progress arises from the interplay of Chiral Perturbation Theory (xPT), the effective theory of QCD at 

low energies [ll-3| and lattice QCD (IQCD) @, in spite of the technical difficulties to perform IQCD simulation for 

(j_) '. light-quark masses close to the physical values. This strategy allows to extract some of the parameters of xPT that 

,.£^ ■ may not be easily accessible in experiments, clarify the role of baryon resonances in the nucleon selfenergy and unravel 

its strangeness content [y, Q- 

A measure of the contribution from explicit chiral symmetry breaking to the nucleon mass is provided by the 
^ \ so-called sigma terms. In particular, the pion-nucleon cr.n-Ar-term is defined as 

m 

00' a^N ^m{N\uu + dd\N) , (1) 

^ _ 

f^ , in the isospin limit to„ = md = m k. A MeV. Using the Hellmann-Feynman (HF) theorem, Gj^j^i can be related to 

^- ■ Mn W^ 

O ■ Q 

m ; a^AT = m-^MN (m) . (2) 

y—^ am 



Additionally, cTt^^^ is the nucleon scalar form factor coming from light quarks at zero four-momentum transfer squared. 

As such, it enters quadratically in the scattering cross section of supersymmetric dark-matter particles with nucleons. 
J, Uncertainties in the determination of sigma-terms, including cr^w, currently represent the largest source of error in 

d , direct dark-matter searches |lll4l3l |. 

Traditionally, the pion-nucleon sigma term has been isolated by extrapolating TriV-scattering data to the (unphysical) 

Cheng-Dashen point [t ~ 2A/^, s ~ u ~ A/^, where s, t and u are the standard Mandelstam variables) [1J| using 

dispersive techniques. The results over the past three decades, Ut^n = 49±8 MeV [151, — 45 MeV [ia|, 56 ±9 MeV |17| . 

64 ± 7 MeV [ll|, 66 ± 6 MeV [Hi, 43 ± 12 MeV [l^]/ depend on the data used as input and on the extrapolation 
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procedure. The lack of consistency among the data sets as well as discrepancies between the parametrizations of the 
experimental data are a sizable source of systematic uncertainties. 

In order to sort the systematic effects out, much effort has been made in the context of baryon xPT (BxPT). At a 
given order in the chiral expansion, B^PT allows to express both the nucleon mass (and ctttat) and the TrA^-scattering 
amplitude in terms of the same unknown low-energy constants (LECs). The available experimental information on 
TrA^-scattering can be used to obtain these LECs. Such a program has encountered a number of difhculties. Unlike 
in the meson sector, in B^PT the power counting (PC) is violated by the presence of Mm as a heavy scale that does 
not vanish in the chiral limit. As a consequence, the loop diagrams do not fulfill the naive chiral order dictated by 
their topology Q ■ The solution to this problem follows from noticing that the genuine non-analytic chiral corrections 
indeed verify the PC, while the breaking pieces are analytic and can be renormalized into the LECs. Different schemes 
have been developed, including non-relativistic heavy-baryon (HB) Q, and the fully covariant infrared (IR) [22| and 
extended-on- mass-shell (EOMS) |23l.[24|. In HBxPT, it was found that the convergence problems in some kinematic 
regions render the fits insensitive to the leading-order contribution to a-j^N- The poor convergence can be traced back 
to the fact that HB limit modifies the analytic structure of the -kN amplitude [23. To overcome the problems of 
HB^PT, the covariant formulations were developed. In the IR approach, loop functions are split into an infrared 
singular part which fulfills the PC and a regular part, containing the PC-breaking terms and higher order ones, which 
is dropped. An important drawback is that the IR scheme introduces unphysical cuts (2^ which can have disruptive 
effects in low-energy phenomenology [2a,l23|- After applying this method to 0{p^)^ Becher and Leutwyler concluded 
that the IR chiral representation of the TrA^-scattering amplitude is a good approximation only in the subthreshold 
region so that no reliable determination of the sigma term could be performed from data in the physical region [28| . 
In the EOMS, the PC is restored by renormalizing the finite number of PC-breaking terms. In this way, the analytic 
structure of the theory is preserved. Two recent EOMS studies of -kN scattering at order p^ [23 and p^ [30| have 
achieved a good description of the data and improved convergence. 

A different complication concerns the treatment of the A(1232) 3/2+ resonance which is only ~ 300 MeV heavier 
than the nucleon and couples strongly to the ttN system. In B^PT, the A(1232) is often treated as a heavy state whose 
influence in the observables is encoded in some of the LECs but, aiming at a more realistic description, it has often 
been taken explicitly into account. In order to include the A(1232) as a degree of freedom one needs to define a suitable 
PC for the new scale A = Ma — Afw |3ll - l33fl . and to treat the so-called consistency problem afflicting interacting 
spin-3/2 fields (see Refs. |34l436| and references therein). The importance of explicitly including the A(1232 ) in By PT 
has been stressed by a recent analysis of the ttN scattering amplitude performed in the EOMS scheme |29l . ISOl l37| . It 
was shown that the inclusion of the A resonance in a covariant framework is essential for a reliable extrapolation to the 
Cheng-Dashen point [33|- The resulting values of CTttAt are in the 40-60 MeV interval, depending on the partial- wave 
analysis used as input and in agreement with with those obtained by dispersive techniques [29J. Although a value of 
cr^jv = 59 ± 7 MeV [221 becomes eventually favored on the grounds of consistency with irN phenomenology [33 , an 
important conclusion of these works is that further efforts are required to understand the possible systematic errors 
in the iiN scattering data. 

Another way towards the determination of the iiN sigma term is provided by IQCD studies. Two different procedures 
have been used. In the first one, the matrix element in Eq. ([1]) is directly obtained and extrapolated to the physical 
values of the quark masses. The second procedure consists of using Eq. ([5]), after a suitable extrapolation of IQCD 
results for Mm down to the chiral limit. The latter has been favored because of the technical difficulties that arise in 
the direct determination of disconnected contributions to ct^at. 

The last decade has witnessed an impressive development of IQCD simulations. Results with two fully dynamical 
light (as light as possible) degenerate fermions [Nf — 2) or with two degenerated light and one heavy (close to the 
physical strange-quark mass) heavy fiavor (N f = 2 + 1) have become standard. Even a direct determination of cr,r7V 
for N j =2-1-1 + 1 (including dynamical c-quarks) has been reported [38| . Baryon xPT provides a natural framework 
to extrapolate lattice data for Mm with heavy quarks down to the physical and chiral limits, provided that the quark 
masses are small enough to warrant its applicability. In the context of HB^PT with a cut-off regularization it was 
already realized that non-analytic terms were important [39 44 ill . The quark-mass dependence of Mm has also been 
investigated with SU{2) IR^PT to 0{p'^) without explicit A 42| and using phenomenological information to constrain 
the input parameters. Baryon xPT also allows to take finite lattice-volume corrections into account, as it was done 
for Mm in Ref. (Hi. A more complete ©(p"*) IR^PT study [Hi included the leading 0{p^) contribution of the A 
resonance with the small-scale expansion and HB approximation. According to this work, the introduction of A(1232) 
as a propagating degree of freedom is not crucial for Mm- This is in contrast with the findings of Ref. [43 made with 
the EOMS scheme up to 0{p^). 

More recently, the M^ dependence of new Nf = 2 IQCD data for Mm has been investigated with HByPT [4g| and 
IRxPT without explicit A degrees of freedom |47l - l49| . The results for aj^M range from 37 to 67 MeV. In the case 
of Ref. [48|, a direct measurement of aT^M |50| was incorporated to the fit, which allowed to increase the precision. 
Furthermore, three new direct determinations of ct^jv have also been performed applying noise reduction techniques 



for a better determination of the disconnected contribution [51| . 

Several collaborations have pursued IQCD simulations of the masses and cr^Ar using Nf — 2+1 configurations |52h62||. 
The extrapolation to the physical point allows to determine g^ at together with other sigma terms and strangeness 
content of baryons. The difficulties encountered in HBxPT [571 . [63fl to accomplish this program were overcome 
applying cut-off regularization schemes |64l.[6a|. using covariant formalisms up to O (p"^ ) fGmGp and 0{p'^) |69l - t73| . or 
complementing HBxPT with an expansion in the inverse number of colors (large- TVp) [22, Hi IZSI • Although SU{3)- 
flavor calculations have reached a considerable degree of maturity, the large number of unknown LECs at O(p^) and 
the size of the current IQCD data set limits, at present, the accuracy attainable in the sigma terms. 

Alternatively, SU{2) B^PT can be used to perform extrapolations of Mn and (TjrW in the light-quark masses with 
the implicit assumption that the influence of the strange quark is embedded in the LECs and that its mass in the 
simulations is close enough to the physical one. The chiral expansion is expected to converge faster than in 5*1/(3) 
BxPT and the different LECs appearing at 0{p'^) can be independently determined using ttN scattering. On the 
other hand, in comparison with the Nf = 2 simulations in which the strange quark is quenched, the extrapolated 
quantities from Nf = 2 + 1 should be closer to those in the physical world. Analyses of Nf = 2 + 1 simulations 
with S'C/(2)-HBxPT ansatzes at 0{p^) and without A(1232) have become standard |63l. ITo. l77|. In particular, with 
HBxPT up to 0{p'^) it was found that aT^N = 84 ± 17 ± 20 MeV with explicit inclusion of the A resonance, and 
CttN = 42 ± 14 ± 9 MeV without it [7g . While the inclusion of the A had little impact on the value of nucleon mass in 
the chiral limit, the central value of the sigma term changed by a factor of 2. It was also pointed out that the lattice 
data exhibited a surprisingly linear dependence on Af^, a feature also shown by other IQCD data [57l |. The importance 
of the A(1232) in extrapolations of IQCD data on Mn has also been recently stressed in an analysis combining B^PT 
and the large- iVc expansion. [Tg ■ Finally, a different strategy was adopted in Ref. [30| , according to which the LECs 
in SU{2) BxPT up to 0{p'^) were determined in simultaneous fits to ttN scattering data and IQCD results. 

Here we present our study of the pion mass dependence of the nucleon mass in covariant SU{2) B^PT up to 0{p'^), 
using the EOMS scheme with explicit inclusion of A(1232) intermediate states. We perform global fits to recent 
determinations of Mjv in IQCD simulations with Nf ~ 2 and Nf = 2 + 1 dynamical quarks, taking into account 
finite lattice-volume corrections. By extrapolating the fits we determine the nucleon mass in the chiral limit and 
the pion-nucleon sigma term, paying attention to the different sources of systematic errors: the extrapolation to the 
continuum of IQCD data with finite lattice spacing, normalization errors, the uncertainties in the LECs fixed in the 
fits and the range of applicability of the chiral expansion. The article is organized as follows. In Sec. |ll] we describe 
the formalism, derive the formula for the nucleon mass and discuss the origin of the different coupling constants and 
LECs that are constrained in the fits. Finite volume corrections and continuum extrapolations are also discussed. 
The fit strategies and the results are presented in Sec. IIIII We conclude and summarize our work in Sec. IIVI Further 
details about the calculation can be found in the Appendices. 

II. NUCLEON MASS IN THE BxPT 

Our aim is to study the pion mass (M^r) dependence of the nucleon mass (Afjv) and obtain the value of the a^^N- 
term by means of the HE theorem. For this we employ the O (p^) covariant SU (2) B^PT with and without explicit 
A-isobar degrees of freedom, AB^PT and /^BxPT. The resulting function AIn (Mt^) depends on several LECs whose 
values we fix by fitting IQCD nucleon mass data for unphysical quark-masses. The required ingredients are established 
in this Section. We derive the perturbative nucleon mass and show the explicit fit formulas together with a discussion 
of IQCD discretization effects. 



To define the nucleon mass in terms of an expansion in the light scales m^ = 2Bm, p and A = A/a — M^ 



N, 



we have to choose a counting scheme^. If A-isobars appear explicitly, the common assumptions are the small-scale 
expansion |3ll . l32J that counts A ~ p ~ m^ and the S — counting [33|, which takes A ~ p^^^ to preserve the hierarchy 
p ~ m^ <C A. As the latter is not the case for most of the IQCD simulations, we adopt the small-scale counting. The 
order n of a self-energy contribution is then defined by 

n = 4L - 2iV^ - A^w - iVA + ^ fc\4. , (3) 

k 

for a graph with L loops, N^^ internal pions, N^ internal nucleons, N^ internal A-isobars and T4 vertices from a £'''') 
Lagrangian. In Fig. [T]we collect all one-particle irreducible diagrams that fulfill, after a suitable renormalization, Eq. 



The constant B is proportional to the chiral quark condensate. 
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Figure 1: One-particle irreducible contributions to the nucleon self-energy up to O (p"*)- Single solid lines denote nucleons, 
double solid lines, A-isobars and dashed lines, pions. Boxes represent the pion-nucleon and contact vertices where the number 
specifies the chiral order. 

([3]) up to n = 4 [©(p^)] and list in App. |2]all relevant B^PT Lagrangians. Among the A-isobar contributions, the 
graphs SArA4a and SArA4h originate from the C\j^^ Lagrangian [23|. It was shown in Ref. [73| for the HB^PT case 
that these couplings are redundant and can be absorbed in the LECs of £^^ and C^lf^ ■ The HBxPT expressions are 
the leading order contributions to covariant B^PT results which implies that these two diagrams start to contribute 
at O (p^) . We do not include them in our O (p^) calculation. Additionally, the ttN scattering analysis [2^ performed 
explicit fits with and without these terms and found strong arguments to support that these redundancies also carry 
over to the covariant case. 

To calculate the remaining diagrams we apply the EOMS renormalization-scheme [2j, l23[which uses the analyticity 
of the power-counting breaking terms to overcome the power-counting problem found in [3|. Explicitly, we calculate 
these diagrams in the dimensional regularization for D = A — 2e dimensions and renormalize terms proportional to 
L = — i — 1 + 7s — ln47r (AfS-scheme). Subsequently, we renormalize the appearing LECs in such a way that 
power-counting breaking terms are canceled. 



A. Nucleon self-energy and the perturbative nucleon mass 



The nucleon physical mass Mjv is defined by the pole position at |i = AIn of its full propagator 

1 



where S (jJ) and Mq are the nucleon self-energy and the (chiral limit) bare mass. In order to define a perturbative 
nucleon mass, we expand E {^) around ^ = Mq: 



E(i^) = j:{Mo) + {^~Mo)-^ 






j:i^) + ly^M,f^^ 



^=Mo 



SW 



and write the propagator as 



1 



^-Mo-^{^) ^_j^j^_ S(^^o) 



Rj^) 1 - S' {Mo) 



1 - S' (Mo) 1 - E' (Mo) 
Equation ^ defines now the nucleon mass by the pole at -^ ~ Mjq 

Mn = Mo + ZY. (A/o) + ZR (Mn) , 

together with its residue 

1 

Z 



1 - S' (Mo) 
Using the B^PT self-energies up to order p'^ of App. [BJ gives: 

Ep4(|i) = E(2) + e(3) (j^) + eW (j^) 

= E(2) + e(3) (Mo) + e(4) (Mo) + {i>'- Mo) [s<'^' (Mo) + Ef^)' (a/q)] + R {^) 
Z = l + E(3)'(A/o)+0(p3) , 



(5) 
(6) 

(7) 

(8) 
(9) 

(10) 

(11) 
(12) 



where the upper indices denote the chiral order. Only the contact term Sc2 = —^Ciirv^ enters in S*^^^ so it does not 
depend on |^. Inserting Eq. (|lip in Eq. ([S]) one gets the nucleon mass up to order p^: 

M^"' {ml) = Mo + Ec2 {ml) + ^m K) + ^naz K) (13) 

+EAr4 {ml) + Et4 K) + SC4 (m^ ) + Sc2 (m2) E^3 (m2 ) 

+ SArA4 K) + Y.C2 K) S^^g (m^ ) + O {p^) , 



where ah loops are evaluated at |i = A/q. The term R{Mn) contributes only at 0{p^). The first line of Eq. (fT3)) 
corresponds to the p^ nucleon mass while the second and third lines are the additional p^ contributions; the notation 
of the different terms matches the one of the diagrams in Fig. [TJ All S^ are obtained from the Lagrangians in App. \K\ 
and are explicitly given in App. [Bj There are 10 low energy constants, namely, /jro, gAo, ci, C2, C3, Haq, Mq, Mao, cia, o.- 
Most of them are constrained by experimental data. More details about their treatment are given below. 



Applying the HE theorem 



to Eq. (fO)) one obtains. 



B. Nucleon mass, CT,rJV-term and fit formula 



<JnN {ml) = rn—AlN (m) - m^^-^A/jv {ml) (14) 



Mlf)(m2) ^ Mo~cAml + ^aml + E\f+^^\ml,Mo,MAo,Uo,9Ao,hAo,c,) , (15) 



d 
drivi 



(4) /_2^ ^ _4cim2+am4+m2-^s(3)+W^^2^^^^^^,^^^^j^^^^^^^;,^^^^ ^ (Ig) 



0■l^^ TO, 



TAf 



with Ci = ci, C2, C3, ciA- The a^^j^ can also be obtained from a direct calculation of the nucleon scalar form factor Eq. 
([T|) at zero four-momentum transfer squared. We have checked that Eq. (J16p can be mapped term by term to such a 
calculation, i.e. that our formulas with full, non-expanded loops fulfill the HE theorem. 

To apply Eqs. p5ll6l) with a p^ accuracy, we cannot identify the physical (or lattice) pion mass Af^ with the lowest 
order m-^ {M^ = to, J = 2Bm) but must take the next order into account. According to the well known expansion [2|: 

AC {ml) =ml + ^-^^mt + —^m^ In ^ + O (/) , (17) 

where l^ (A^) is a renormalized scale-dependent LEG coming from the meson xPT Lagrangian. Therefore 

Mj^'> {Ml) = Mo - c^AMl + ^mi^ + ^j^Mt In ^ + lfj+^^^ {M^ Mo, A/ao, U, gA, Ha, q) + O {p') , (18) 



a. 



a {Ml) = -4crA42 + aAf;J-crA^3^(M„2)A/;J + ^Af;Jlng 

+Mlg^^'i!l^l'^ {Ml, Mo, Mao, U,gA, Ha, q) + O (/) (19) 

16 

with a = a + ci — ll{Ml) . (20) 

J-n 

Equation ^ is our final formula for O (p^) BxPT fits to IQCD data. The effect of Eq. ^ is an additional O {p'^) 
term proportional to ci and a redefinition of a — )■ a which will be a fit parameter. Furthermore, we adopt the physical 
values of /^ = 92.4 MeV and gA = 1.267 instead of the chiral limit ones and set the renormalization scale to A = Afo. 
The differences between the chiral limit and physical values are of order p^ so that they start to contribute at O {p^) . 

In the case of a\li {Ml) we cannot absorb all terms proportional to /g (A/g ) in the LECs and shall need a numerical 
value for it. From the latest estimate of Is (A/^) = InA^/A/^ at the physical point Is (139 MeV) = 3.2 (8) 0, H^l one 
has: 

H (A^) - -^ (^3 (M.) + m §) ^ -^ (3.2 (8) + m %^) , (21) 

where we set M^(^pi^yg) = 139 MeV. 



Theory 


C2 [GeV-^] 


C3 [GeV"''] ciA [GeV-^] /ia Mao [MeV] 


4^-BxPT 
A-BxPT 


3.9 ±0.4 
1 l+°-2 


-6.7 ±0.4 
-3.0l!^;^ -0.90 ±40 2.87 1170 ±30 



Table I: Values ol the LECs appearing in the p'' nucleon mass. For the LECs fno and qao we take their physical values /^ — 92.4 
MeV and qa = 1.267. 

C. Low-energy constants, finite volume and lattice spacing effects 

After fixing f-^Q and qao, we discuss our treatment of the remaining eight LECs, ci, C2, C3, Mq, a, Mao, /iaoj 
and ciA- Generally, our fits depend very mildly on variations in C2, C3, Mao a-^d ^^A- Furthermore, we observe that 
changes in cia are compensated by changes in a. Our strategy is, therefore, to fit Mq, c\ and S while keeping 02, cr, 
Mao, ^a and cia fixed. The nucleon-related LECs C2 and C3 are taken from the TrA^-scattering analysis of Ref. [37|, 
performed with the same BxPT framework employed here. More specifically, we take as central values the average of 
the results of fits to the phase shifts from the Karlsruhe-Helsinki group (KA85) and the George Washington University 
group (WI08), accepting errors defined by their uncertainties and also by the result of the fit to Matsinos phase shifts 
(EM06) (see Tables 1 and 2 of Ref. [13) ^. The specific figures for both the ^-BxPT and A-BxPT cases are given in 
Table H 

In order to fix the A-related LECs, Mao, cia and Ha, we consider the pion-mass dependence of the A-isobar mass. 
Up to O (p3) it reads [Hi 

Mf ^ (M,) = Mao - 4ciaA/^ ± ^an3 (M,; Ha, U, Mn, A/a) ± Saa3 (M,; Ha, U, Ma) , (22) 

where the loop contributions Sans and SaA3 stand for diagrams like S^ra and SjvA3 in Fig [1] but with external 
nucleon lines replaced by A(1232) ones. The explicit expressions are given in App. |B] Since differences are of higher 
order, we are allowed to take phenomenological values for the LECs in these loops, instead of the leading-order ones. 
In this way, one uses the phenomenological value of the A-isobar decay width TA-^N^T = ^ 2ImI]AJV3 = 115 MeV to 
fix Ha = 2.87. Furthermore, we adopt Ha = IffA obtained in the large-A^c Hmit. Finally, we use IQCD data for the 
A(1232) [5^; [53, [81| mass to determine the remaining two LECs Mao and cia. As the available lattice results are 
rather scattered, we do not perform a rigorous fit to them but, instead, adopt the conservative attitude of setting a 
band that englobes all the IQCD points with their errorbars (see Fig. [2]). The central values for the parameters result 
from the average of those defining the band's boundaries and are listed in Table H] 

We now turn to two discretization artifacts: finite volume (FV) and finite spacing effects, appearing in IQCD 
studies, as a consequence of the finite grid with volume L^ and spacing a in which simulations are performed. 

All loop graphs of Fig. [T] are subject to FV corrections. We calculate them in App. IB 31 applying the standard 
techniques of Ref. ^3|- The FV corrections to Sjvs and Et4 are equivalent to those in Ref. |4J|. In addition, we 
correct the combination Ejv4 + Ec2S^3, the A-isobar graphs Satas, ^NAi and Sc25^jvA3 which contribute at order p^ 
in the continuum [2j,[2J| in the EOMS renormalization scheme. Reference [43 employs IR, for which the combination 
of S7V4 ± Ec2S'^3 appears only at order p^. Our FV corrections are therefore: 

Ep4 {Ml L) = Sats {Ml L) ± SArA3 {Ml L) 

±Sw4 {MlL) + S]^A4 {MlL) + ST4 {Ml L) 

±Sc2 {Ml) S^3 {Ml L) ± Sc2 {Ml) S^^3 {mI L) . (23) 

All these terms are given in App. IB 31 In Fig. [3] we test our FV correction against IQCD data with approximately 
the same pion mass but different L. We found four points from the QCDSF collaboration [43, four points from the 
NPLQCD collaboration [6l| and two points from the ETMC collaboration [S^l at pion masses approximately of 265, 
300, 390 and 440 MeV, respectively. Reasonable values of the LECs A//o = 890 MeV and ci = cia = -0.9 GeV'^ have 
been chosen for this exercise. We observe that our FV corrections describe very well the L dependence for lattice 
sizes larger than ^ 2.2 fm and that they have a size of up to 45 MeV. In our fits we shall include only data points 
with LM^ > 3.8 for all of which L > 2.2 fm. 

In general, we will use IQCD data that are not extrapolated to the continuum limit a — > 0. Originally, discretized 
QCD actions break chiral symmetry even in the chiral limit by terms proportional to a |83l - l85j but modern lattice 



Further justification for this choice is given in the Results Section. 
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Figure 2: Pion mass dependence of the A-isobar mass. 
Green squares are from 152 . \53^ . black right-triangles 
are quenched data from [8ll | and red diamonds are un- 
quenched data from |8ll |. The blue circle is the physical 
point. The band defines the uncertainty range adopted 
(see the text) while the blue line is the preferred result. 



Figure 3: Finite volume corrections AAfjv — A4n{L) — 
Mn{L -^ cxo) as a function of the lattice size for pion 
masses of 265 MeV, 300 MeV, 390 MeV and 440 MeV. 
Lattice data from Refs. |4a | (triangles), |6ll | (red dia- 
monds) and [83 | (squares) with approximately the same 
pion masses are also displayed. We normalized each curve 
to the point with the largest volume and shifted them by 
multiples of 50 MeV to avoid overlaps. At L = 4.0 fm 
AMjv ~ for all curves. 



calculations use O (a) improved actions for which discretization effects in baryon masses start at order a^ . However, 
there exists a whole variety of IQCD-actions, each with its own discretization effects. For the specific Symanzik 
IQCD action an effective field theory investigation has been performed in Ref. J8^ on a HBxPT basis but a general 
approach, similar to the treatment of FY corrections, does not exist. Therefore, we parametrize this effect for each 
action individually by writing the nucleon mass in an a -expansion to the lowest order as 



M 



N 



M, 



a=0 



+ Caa^ + O (a' 



3 2 2^ 

a m^ 



(24) 



with an action-specific constant Ca- By using the ETMC points at A/^ = 260 and 262 MeV, and QCDSF points at 
toMtt = 0.658 and 0.660 |48l. |83| we can roughly estimate the size of this effect. By taking the linear a^-extrapolation 
of Eq. (PH) we obtain Cetmc = 0.17 GeV'^ and cqcdsf = 0.33 GeV^, which correspond to nucleon-mass shifts 
of 10 — 50 MeV. We obtain that lattice spacing corrections can have similar sizes as the FV ones. Therefore, we 
incorporate this effect in specific fits by including the CaO^ term in the x^ for each collaboration/action reporting 
results for different values of a. 



III. RESULTS 



We study the pion mass dependence of the nucleon mass by using the covariant B^FT expression of Eqs. (|13p and 
([T5)l , which is accurate up to the chiral order p^ and includes explicit A-isobar degrees of freedom. We perform global 
fits to IQCD ensembles for Nf = 2 and Nf = 2-1-1 numbers of fiavors. Generally, IQCD uses a discretized QCD-action 
to simulate the quark- gluon interaction in a finite box of size L^ xT with finite spatial and time spacings of a and a*. 
The nucleon mass data are given in terms of the dimensionless quantities aM^^ and aM^ with uncertainties in a, aM^^ 
and aMpf. An actual value of a sets the overall scale to convert the IQCD data into physical units. No universal scale- 
setting method exists and different collaborations use different approaches. Furthermore, the statistical uncertainty in 
a turns into a normalization uncertainty in M^ for data points belonging to the same a-set. It is therefore preferable 
to fit the [aM-^^aMiq) data directly whenever this is possible or, otherwise, to include these correlated uncertainties 
in the fit. As explained below, we are able to perform the former in the case of the Nj = 2 ensembles and rely on the 
latter for the Nf = 2-1-1 ones. We also include FV corrections and lattice spacing effects as described in the previous 
section. We fit the LECs Mq, ci and a while keeping C2 ,03, cia, hA and A^ao fixed to the values listed in Tab. H] 
Afterwards, we quantify the effect of varying the fixed LECs within their ranges. The fit uncertainties are determined 
at a 68% confidence level. 
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For Nf = 2 we include data from the BGR [sljETMC [si^ainz [H and QCDSF jii collaborations, and for 
Nf =2 + 1 from the BMW Q, HSC [Hi, LHPC [5|, MILc|8i|, NPLQCD [Ml, PACS @ and RBCUK-QCD ^ 
collaborations. In both cases we extract the LECs and obtain the aT^N value by using the HP theorem. 

A. Nucleon mass up to order O (p*): fits to Nf = 2 lattice QCD data 

We use Eg. (fTSl) to fit the IQCD data for the Nf = 2 ensembles of the BGR, ETMC, Mainz and QCDSF 
collaborations [481 I82I . 1871 . I88| . The IQCD data are given in terms of the dimensionless products aM^r and aM^ where 
the scale is fixed in different ways: with the experimental fl^ mass in Ref. [83| and with HBxPT or IR-xPT chiral 
extrapolations of A/jv in Refs. [4^ i82|. The available information for these data sets is such that we can perform our 
own scale setting. By doing this we compensate for the different scales of the various sets and avoid to manipulate 
them with two different B^PT versions. 

Explicitly, we fit the IQCD data in terms of (roA/7r,roMAr) by using the Sommer-scale rg [9l|| and the ratios ro/a 
in the chiral limit, as reported by each Collaboration. The uncertainties in uMt^, aM^ and r^/a are assumed to be 
uncorrelated. The value of Tq is a priori unknown and we determine it recursively inside the fit. This is the same 
strategy used in Ref. [i^, now employed to analyze Nf = 2 data globally. The x^ function that we minimize is 



x^ = E 



{mi) 



a4'^ ( Af 2) + sj;') ( A/^, L ) + caa' - d, [ M^, L 



(25) 



with A/];^=roA4"^ , A4 = {roA'Uf , S^^ = roS^;^ , (26) 

where di I M^,L] are the IQCD data points with uncertainties Ci, each of them generated in a lattice of size L and 

spacing a. The continuum expressions Mj^ (M^) and the finite volume corrections E^"^ {^M^,L) for the chiral-order 

n are listed in App. IB 41 As discussed above, the terms CaO,^ = r^Ca [a/ro) parametrize discretization effects, with Ca 
being common constants for points obtained by the same IQCD Collaboration/action. The Sommer-scale is calculated 
in each minimization step recursively using the constrain imposed by the experimental value of the nucleon mass at 
the physical point: 

MN(^phys) 

The explicit fit parameters in Eq. (pS)) are Mq, ci, a and two Ca constants, one for the ETMC Collaboration and one 
for both Mainz and QCDSF which employ the same action. The single data point of BGR does not allow to perform 
any lattice spacing correction. As the term CaO/^ does not stand on the same firm ground, from the perspective of 
effective field theory, as the rest of our mass formula, we perform fit with and without it and treat the differences as 
systematic errors. We restrict the data sets by imposing the following conditions: roAf,r < 1-11, Af^L > 3.8, which 
englobe points of A/jr < (429,476) MeV for Sommer-scale values in the range ro = (0.51,0.46) fm. We then consider 
the following data sets 



• 



BGR [83 : A Sommer-scale of tq = 0.48 fm is assumed and three data points are provided, only one below 
toMn = 1.11. 



• ETMC [82|: Eleven data points are provided in the form {aMT^,aMjf)] for each setting a value of r^/a is 
computed. After converting {aMTr,aMpf) into {toMtt, tqMm) we find that seven data points fulfill our conditions 
and enter the fit. 

• Mainz [8^: Eleven data points are provided in the form (aAf^, aMjq). The lattice spacings as well as the ratios 
ro/a are determined by the Vt~ mass [92,[93|. We convert (aMT^.aM^) to (roAf7r,''o-^Af) and six data points 
enter the fit. 

• QCDSF [43: This work provides 27 data points, directly in terms of (roAf^, roA-fAr), but only two of them 
fulfill our restrictions. In addition, there is a single data point for the a^^N obtained by direct determination at 
AU -- 285 MeV [13 . 

We study the following variations of the fits: 

1. Mis!{M.,r) to order p^, p^ and p* in the chiral expansion 







excluding 


a-^N (285 MeV) 






including 


a^N (285 MeV) 






Mo [MeV] 


ci [GeV-'l Q 


[GeV-3] 


Lf '■o M 


a. [MeV] 


Mo [MeV] 


ci [GeV"'] a 


]GeV-3] 


dof '^0 IH 


a^ [MeV] 


p' 


906(11) 


-0.43(2) 


- 


2.1 0.509 


34(2) 


913 (6) 


-0.33(1) 


- 


6.3 0.539 


26(1) 


p' 


880(13) 


-0.93(3) 


- 


1.9 0.480 


53(2) 


892 (6) 


-0.78(1) 


- 


8.5 0.527 


41(1) 


pi 


863 (16) 


-1.19(4) 


- 


2.1 0.456 


68(3) 


878 (5) 


-1.00(1) 


- 


9.5 0.517 


52(1) 


p' 


866 (40) 


-1.18(14) 


23(3) 


2.5 0.470 


62(13) 


888 (9) 


-0.91 (4) 


38(2) 


2.9 0.507 


41(3) 


pi 


893 (29) 


-0.77(9) 


35(2) 


2.4 0.494 


38(10) 


890 (7) 


-0.80(1) 


33(2) 


2.5 0.489 


41(2) 



Table II: Results for BxPT fits to Nf = 2 nucleon mass data from Refs. [43, [83, 113, JSQ] • The 'A' index denotes the inclusion of 
explicit A- isobar (ABxPT), while its omission corresponds to /^BxPT; FV corrections are included but finite-spacing effects 
are excluded. The left-panel results come from a fit of solely nucleon-mass data while in the right panel the a-^N point at 
M^ — 285 MeV of Ref [501 was also taken into account. 
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Figure 4: Fits to the Nj — 2 nucleon mass data of Refs. [4a. |82. ISTI. ISq . Filled (open) symbols are for data points included in 
(excluded from) the fits. The left (right) picture shows fits without (with) explicit A-isobar. The fit including the cr^jv (285 MeV) 
of Ref. |50[ is given by the blue solid line while the plain nucleon mass fit is given by the green dashed one. The dark blue 
and light green shaded regions represent the corresponding statistical uncertainties. The IQCD data are scaled by ro and FV 
corrected according to the simultaneous fit. Hence, the green dashed line does not correspond to the shown data points. 



2. without (^BxPT) and with (ABxPT) A-isobar 

3. including and excluding the single direct cr ttA? measurement of Ref. [50| 

4. without and with lattice spacing corrections [caO? term) 

5. variations of the input LECs according to the errors quoted in Table [J 

Finite volume corrections are always included. 

The output of our fits for cases 1-3, with the LECs fixed to the values in Table [J and without lattice-spacing 
corrections are presented in Table [lT| and Fig. [31 Bear in mind that changes in the fit conditions 1 and 2 yield different 
To (see Table [nl) so IQCD data are scaled differently. From Table [Til we observe that the inclusion of O(p^) does 
not lead to a better description of present nucleon mass data than the 0{p^) one. However, for fits including the 
(Jtjn (285) point, a good \^ /dof emerges only at 0{p^). In this situation, ABxPT gives a slightly better \^ /dof 
than i^BxPT but both approaches give the same (TttW value. The overall rather high x^ /dof is caused by two points 
from the Mainz Collaboration. By excluding them we obtain -)^ /dof ~ 1.6 but the results change only within the 
quoted uncertainties. The FV corrections shift the data points by (—6) — (—50) MeV. In contrast to the z^BxPT case, 
the ABxPT p^-results are not significantly altered by the inclusion of Ot^j^ (285) in the fits and exhibit a softer Mjr 
dependence. This might be interpreted as an indication that the theory with explicit A(1232) is more realistic. 

Figure El shows the relative contributions, \j? /'p^\ and jp^/p'^j, of different chiral orders to the nucleon mass for fits 
including o-ttTV (285). One observes that the O(p^) term has a relatively small contribution over a large Af^ range. The 
same is true for the :;2^BxPT 0(p'^) term. In the ABxPT case, however, the relative impact of the Oij?) contribution 
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Figure 5: The B^PT results for MN{M.n) decomposed 
into their chiral-order relative contributions |p"^/p^| and 
|p*/p^|. The blue solid line denotes |p'^/p^| and the 
purple-dashed line, b^/p^l, both for ABxPT. The red 
dashed-dotted and orange dotted are the [p'^/p^| and 
|p*/p^| results for If^'x^'V. The shaded region is excluded 
from the fit. 



Figure 6: Pion-mass dependence of the cr^Ar-term. The 
blue solid line and the green dashed lines stand for the 
AB^PT and .^ByPT fits respectively, both including 
a-^fq (285) of Ref [501 (red diamond). The dark blue and 
light green shaded areas represent the corresponding un- 
certainties. The red square is our final result at the phys- 
ical point. 



pi 


A/o [MeV] ci [GeV-i] a [GeV-=*] ce [GeV-^*] cmq [GeV-^'] ^ n, [fm] a, [MeV] 


excluding cr,riv (285 MeV) 
including a^iv (285 MeV) 


894(28) -0.76(10) 36(5) -0.06(7) -0.05(13) 2.8 0.501 37(10) 
892(21) -0.79(2) 34(3) -0.08(6) -0.08(12) 2.8 0.499 40(3) 



Table III: Results for p^'-ABxPT fits to iV/ = 2 nucleon mass data from Refs. [4|,|82 
by the ceo,^ and cmqo^ terms for the ETMC and Mainz/QCDSF data respectively. 



with lattice spacing effects accounted 



steadily rises, becoming more than 80% of the p^ one at Af^ > 450 MeV. From this we deduce that A/.^- ~ 450 MeV 
is at the upper border of the AB^PT applicability. We have also performed fits with relaxed conditions LM^^ > 3.5 
and toMtt < 1.00 which, however, yield equivalent results to those already presented in Table HIl The present data 
does not allow to go below roM^r < 1.00. 

In Table IIIII we summarize our results including finite-lattice spacing corrections in the fit, namely the CEd^ and 
cmqh^ terms for ETMC and Mainz/QCDSF respectively. We obtain corrections of (-1-6) — (+20) MeV, which have 
an opposite sign with respect to the FV corrections. By comparing to Table HIl we notice that all changes are within 
the already given uncertainties. A noticeable qualitative effect is that changes in the Sommer-scale counterbalance 
finite-lattice spacing corrections so that the results are close to the former ones. A more elaborated EFT background 
is required to calculate and interpret finite-lattice spacing corrections more reliably. 

We have tested the fits for variations of C2, C3 within the errors given in TableHl In all cases the results are compatible 
within uncertainties with those of Table |IT1 We conclude that the p'^ B^PT fits are not able to constrain these LECs 
effectively. Furthermore, by varying cia we find it to be correlated with a. The inclusion of cia as a free parameter 
does not produce reliable fits unless the ct^at (285) is taken into account. The fit is driven to unreasonable high cia 
values with rather large a values. However, in fits including the aj^N (285) point we recover cia = —0.87 (30) GeV~^ 
together with results compatible with those in Table [TTl A scan over a range of cia shows that reasonable fits can 
only be obtained in the interval cia — —0.6 ^ —1.3 GeV~^. The better nucleon mass data do not allow to constrain 
ciA more accurately than the present A-isobar mass IQCD data do (see Fig. [5] and the corresponding discussion). 

As a final o-TrAf-value for the Nf ~ 2 IQCD fits we quote 

a^jv = 41 (3) (1) MeV , 

which corresponds to our p'^ /^ and AB^PT fits including a^^N (285) and FV corrections. As for the statistical 
uncertainty we take larger one of ^^B^PT. The difference with the result including finite spacing effects is taken 
as a systematic uncertainty. Note that by excluding the point ct^^t (285) we get a larger cTttN = 52 (13) (11) MeV. 
This value is obtained as the average of the two AB^PT and ^^B^PT results, with the statistical uncertainty being 
the largest of the two errors and the second, systematic uncertainty stems from the difference between AB^PT and 
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:^BxPT central values. 

Figure [S] summarizes our results for the pion mass dependence of the crTrAf-term. The results for the AB^PT and 
/^BxPT fits are compatible within errors but exhibit a different M^r dependence. 

For our final values of the LECs Mq, c\ and a we quote those of the p^-AB^PT fit of Tab. [ll] including (t(285). 
In particular, in the present work we set the Sommer-scale to tq = 0.493(23) fm, which is the average of all our p^ 
results and where the uncertainty is chosen such as to cover all our p'* fits. 

B. Nucleon mass up to order O (p*): fits to Nj =2 + 1 lattice QCD data 

We use our ByPT nucleon mass formula of Eq. p^ to fit the IQCD data for the N f = 2+1 ensembles of 
Refs. |56l . ISSl . ISOl . Wn . Wa . [S^, ISOl- Most of the data are provided in terms of (aMT^^aMN), together with the individual 
lattice spacings a and the statistical uncertainties for all the three quantities. Unlike the Nf = 2 case, the available 
information does not allow us to perform our own scale setting. Therefore, we treat the a-uncertainties as correlated 
normalization errors for all Mn points from the same set. Our treatment of normalization uncertainties follows from 
Ref. j94l |. We perform three types of fits: 1) neglecting correlated normalization errors, 2) including the normalization 
error in scale factors fi, 3) including the normalization uncertainty in a correlation matrix V. For the case 3) we also 
consider lattice spacing effects. The x^ functions for type 2 and 3 fits read 



xl 



xl 



E 



MP {Ml) 



.(«) 

^N 



{MlL)-fA{MlL) 



ftCT'l 



A^ V~'A 



with 



A,, 



M 



(") 

TV 



Mi 



/.-I 



"^f. 



- Cia^ + S^ 



M,^,U 



[M'r,U 



(28) 



(29) 



where M^"^ (Af^) 



and S]^ (Af^,L) are the BxPT continuum and finite volume expressions given in App. 



di (M^, L) are the IQCD data, each point for a given lattice size L and spacing a. We denote the statistical uncertainty 
for Afjv coming from aAfjv as Ui and the normalization uncertainty coming from a as c/. . Case 1) is recovered from 



Eq. 



by taking all /,; = 1 and replacing ai — >■ af 



(T? corresponding to the assumption that Ui and tr/. are 



uncorrelated errors. In case 2) the fi are additional fit parameters; uf and cr?. are treated separately. In case 
3) af and cr? are incorporated in the correlation matrix V . The BMW collaboration [53| does not provide enough 
information to disentangle the uncertainties from aA/^v and a so that we always include this data set with uncorrelated 
uncertainties. 

We include in the fits data points with M.^L > 3.8 and M^^ ^ 415 MeV. The LECs C2, C3 and cia are fixed to 
the values given in Table HI There are two points with M^r ^ 390 MeV from Refs. [5a, |6l| with very small reported 
CTi and slightly smaller Afjv values compared to the neighboring points (see Fig. [7]). The inclusion of these points 
shifts the results to lower masses, yielding a slightly worse x^/do/. Although these points were obtained by different 
NPLQCD and HSC Collaborations, they are not entirely independent because NPLQCD uses the scale of the HSC 
Collaboration, which actually expresses some concern about the quality of their lattice-spacing determination. In 
view of the situation, we exclude these two points from our main results but consider their infiuence in the systematic 
uncertainties. 

In Table HV] we display our results for the fit types 1) and 3). The results obtained with option 2) are similar to 
those obtained with 3) so we do not show them. The consideration of normalization uncertainties slightly enhances 
the x^/do/ but causes a noticeable reduction of ci and aT^N- The quality of the fits in terms of x^ /dof is essentially 
the same for p^ and p'* fits. As in the Nf = 2 case, we expect the advantage of the p"^ formula to be tangible as soon 
as direct CTrTV-data for low pion masses become available for Nf = 2 + 1. 

The left panel of Fig. [7] shows the pion mass dependence of our 0{p'^) nucleon mass results for both /^BxPT and 
ABxPT. There is a large overlap of the corresponding error bands, which are addressed below in more detail. By 
decomposing the fits into their chiral-order relative contributions (right panel of Fig. [7]) , we observe a similar situation 
as in the Nf = 2 case. Namely, the 0{p^) relative contributions are small over a large range of Af^r but the 0{p^) in 
ABxPT increases, making the applicability of our perturbative expression questionable for high Af^ values. We have 
checked that a fit constrained to Af^r < 360 MeV produces results compatible with those of the Af^r < 415 MeV fit but 
with larger uncertainties. 

The results of the fits taking into account lattice spacing effects are given in Table [V] These are considered for 
data sets with enough points with the same L and different a values. Explicitly, we introduced two terms cmcl^ and 
cro^ for the MILC and RBCUK Collaborations, respectively. In the case of BMW, we assume that lattice spacing 
uncertainties are included in the errorbars. We find nucleon mass shifts of (—7) — (—46) MeV, which are small but 
comparable in size with the FV corrections. With this correction, the x^ /dof is slightly better and ct^at gets smaller 
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Figure 7: Combined fits to IQCD data of the iV/ = 2 + 1 ensembles [53, S |5i, |6l|, 162, IH, |90| ■ Left: Fits to nucleon mass data 
up to M^ = 415 MeV. The blue solid (green dashed) line shows the C'(p*) AB^PT (^^BxPT) fit of type 3). The red dotted 
line is also for ©(p*) AB^PT but including the two points of M-^ ~ 390 MeV, excluded from the main fits as discussed in the 
text. Filled (open) symbols represent point included in (excluded from) the fits. Right: Decomposition of the fit results in 
their chiral order contributions. The blue solid line corresponds to the |p''/p^j ratio and the purple-dashed one to |p'*/p''|, both 
for ABxPT. The red dashed-dotted and orange-dotted are the |p'*/p^| and |p''/p'*| results obtained with :^BxPT. 



LM^ > 3.8 
M^ < 415 MeV 


Mo [MeV] 


^BxPT 
ci [GeV-'l a [GeV^] 


dof 


a^ [MeV] 


A/o ]MeV] 


ABxPT 
ci ]GeV~'] a ]GeV-='] 


dof 


a^ ]MeV] 


P' 


904 (2) 
883 (2) 
870 (3) 


-0.47(1) 
-0.90(1) 
-1.15(3) 


24(2) 


3.1 
1.3 
1.3 


36(1) 
51(1) 
58(3) 


870 (2) 
883 (3) 


-1.10(1) 
-0.89 (3) 


26(2) 


1.2 
1.4 


60(1) 
49(2) 


no correl. p^ 
no correl. (390) p* 


865 (5) 
863 (5) 


-1.22(5) 
-1.25(5) 


19(4) 
15(4) 


1.0 
1.4 


63(3) 
64(3) 


878 (4) 
876 (4) 


-0.96 (4) 
-0.99 (4) 


20(4) 
15(3) 


1.1 
1.6 


54(3) 
56(3) 



Table IV: Combined fits to the Nf = 2 + 1 IQCD ensembles [53, \M, \M, M, M, M, W for Pion masses M„ < 415 MeV. The 
LECs C2, C3 and cia are set to the central values given in Table IT] FV effects are included while a^ effects are excluded. The 
Icist two rows correspond to fits of type 1) neglecting correlated normalization errors. The fit of the last row takes into account 
the two points of Refs. [53, l6l|] with M^r ~ 390 MeV, excluded from the main fits as discussed in the text. 

by several MeV. The uncertainties for the constants cmi^ and cro^ are now slightly smaller than in the 2 flavor case 
although all values of Tables Hill and IVl agree within the individual errors. 

We tested our results for changes by varying C2, C3 and cia within the erros quoted in Table [H All changes are 
within the above quoted uncertainties. In particular, changes in cia are compensated by changes in a and reasonable 
results are only obtained for the range of cia = (—0.5) — (—1.3) GeV^^ estimated above. 

As a final value for o-ttW hi the Nf = 2 + 1 case we give 

a^N = 52 (3) (8) MeV , 

obtained in the following way. The central value is the average of the four C(p'*) AB^PT and i^BxPT results without 
(Table lIVp and with (Table |V]) lattice spacing corrections, all including correlated normalization uncertainties. The 
first error correspond to the largest statistical uncertainty of the values under consideration and the second is the 



Mo ]MeV] ci ]GeV-'] a [GeV-^] cm [GeV-^] cr ]GeV" 



X 
dof 



]MeV] 



^BxPT 
ABxPT 



873 (4) 
887 (3) 



-1.10(5) 
-0.84 (4) 



27(3) 
29(3) 



0.18(8) 
0.21 (8) 



0.03 (2) 
0.04 (2) 



1.2 
1.2 



55(3) 
44(3) 



Table V: Combined fits to the Nf =2-1-1 IQCD ensembles [Sy, |62, |8y] including CaO^ corrections for the MILC (cm) and 
RBCUK (cr) Collaborations. The LECs C2, C3 and cia are set to the central values in Table U 
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Figure 8: Pion mass dependence of Mat and cr,riv given by different ©(p*) BxPT fits to Nj =2 + 1 data. The blue solid and 
green dashed lines stand for AB^PT and ^BxPT. The red dotted line is the ABxPT solution with data points only up to 360 
MeV. The black dashed-dotted line does not take correlated normalization uncertainties into account. The blue circle is the 
phenomenological nucleon mass and the red square is our ct^jv result at the physical point. 
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Figure 9: Pion mass dependence of the nucleon mass. The blue solid line and blue squares correspond to our fits to Nf = 2+1 
IQCD data. The red-dashed line and red-triangles correspond to our fits to Nf — 2 IQCD data including the (j,r]v(285) point. 
The errorbands for our fit results have been removed for the sake of clarity. 



largest difference among them. 

Further conclusions can be extracted from Fig. |5] where the pion mass dependence of M^ and a^rN is shown for 
various p'^ fit strategies. We can see that the small slope variations in Mp{{Mt^) (left plot) translate into changes in 
cr^iv of less than 10 MeV at the physical point (right plot). One also notices that the uncertainties of the individual 
IQCD data points (see Fig. |V| tend to be larger than these variations. We do not expect that with more low-A/.^- 
nucleon mass data points one would be able to reduce the Gt^m uncertainty much further, although simulations using 
one lattice action and different lattice spacings would be very important for a systematic treatment of discretization 
uncertainties. On the other hand, Nf = 2 + 1 direct measurements of a-^N at low M^r < 300 MeV would probably 
lead to better constrained fits as it happens for Nf = 2, reducing uncertainties significantly. 

A remarkable outcome of our analyses is the slight disagreement between the determinations of a-^N using either 
Nf = 2 or 2 + 1 data. Indeed, only one out of the four different Nf = 2 + 1 fits used to obtain our final result is 
consistent with one of the Nf = 2 fits. One may argue that the physical content in terms of the strange quark is 
different in the two types of simulations. In 5[/(2)-BxPT (in the present study or in the description of nN scattering 
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data), the contribution of the strange quarks is integrated out into the LECs. For the 2-flavor ensembles, where the 
strange quark is quenched, one assumes that the effect of the strangeness in the nucleonic observables is neghgible, 
whereas one includes its dynamical contributions in the Nf = 2 + 1 case. Intuitively, simulations with 2 + 1 flavors 
should be closer to the real world. 

However, we find it difficult to address this question with the current IQCD results, and to really tell if the observed 
difference could be produced, instead, by some systematic effect not properly accounted for by our analysis. In Fig. ^ 
we plot together the data for 2 and 2 + 1 flavors, as well as several fits obtained in B^PT at p'^. The Nf = 2 + 1 
data are more homogeneously distributed over the low M^r region, where one expects the chiral extrapolators to be 
more reliable. For Nf = 2 there is also a direct calculation of (t-^n in the chiral regime, which influences the fits 
considerably. Moreover, the scale-setting procedure is different in both cases. The physical nucleon mass was used to 
fix the scale in the 2-flavor case, thanks to the fact that the collaborations reported chiral-limit values for the ro/a 
ratios, whereas in the 2 + 1 case we had to rely on independent determinations of the lattice scale and account for 
normalization uncertainties by incorporating correlations in the x^.* In conclusion, we think our analysis exploits 
the considerable size of the current data set on M^ in a way that it is possible to become sensitive to unexpected 
systematic effects. However, more IQCD data will be required to settle this issue and interpret possible discrepancies 
of this type. 

IV. SUMMARY AND CONCLUSION 

We have studied the nucleon mass and the o-TrAr-term in the SU (2) covariant B^PT up to the chiral order p'^. We 
have performed fits using BxPT with and without explicit A-isobar degrees of freedom to combined IQCD data from 
various Collaborations for Nf = 2 and Nf =2 + 1 numbers of flavors. Special attention has been payed to the different 
sources of uncertainties in the input data. This study is the first application of the p^ SU (2) covariant BxPT with 
the EOMS renormalization scheme and consistent treatment of the A-isobar to IQCD data. We have included finite 
volume corrections and also discussed finite spacing effects. In the Nf = 2 case we were able to set the IQCD data 
normalization via the Sommer-scale rg and also performed simultaneous fits to nucleon mass data and one available 
low M^ (Tttjv data point. In the Nf = 2 + 1 case we took into account correlated normalization uncertainties for points 
belonging to the same data set. In the following we summarize our findings. 

• Our formula for the nucleon mass depends on several low energy constants, some of which have been fitted to the 
IQCD data. Explicitly, the LECs are Mq, ci, C2, C3, cia, M^q, gA, fn, h-A and a; the latter is a linear combination 
of several couplings that appear in the chiral Lagrangian at 0{p*). We adopted the phenomenological values 
for gA, fix and Ha- Our fits are insensitive to the chosen values of C2, C3, cia and Mao so that we are 
not able to constrain C2 and C3 and fixed them to phenomenological values extracted in from TriV-scattering. 
Furthermore, we observe that cia and a are correlated, which hinders a better determination of cia than the 
range cia = (—0.5) — (—1-3) GeV^^ based on rather scarce IQCD data for the A(1232) mass. The LECs A^o, ci 
and a are better determined, and their values are listed in Tables [III and IIVI for the Nf = 2 and Nf = 2 + 1 fits. 
For the Nf = 2 ensembles we were able to extract the Sommer-scale, flnding tq = 0.493 (23) fm. By performing 
flts to nucleon mass data alone as well as including a cr-^M IQCD data point at Mjr = 285 MeV from the QCDSF 
Collaboration we have obtained that the inclusion of the p"* order improves the quality of the simultaneous fits. 

• For both Nf = 2 and 2 + 1 ensembles we have investigated the effects coming from finite lattice spacings a and 
volumes employed in IQCD. We parametrized lattice-spacing effects by linear a^ terms. For FV we applied the 
standard BxPT FV corrections. We have obtained that both effects yield comparable numerical corrections to 
the nucleon mass. However, we also found that the simple parametrization of the finite lattice spacing effects 
does not allow to disentangle it in a quantitative manner from other effects. Fit results with and without finite 
a^-effects are compatible within the statistical uncertainty. In contrast to the a^-effects, the FV corrections are 
very better under control by the established BxPT techniques for the presently available IQCD volumes. 

• We have extracted the CTTrAr-term for the Nf = 2 and Nf = 2 + 1 IQCD ensembles obtaining a-^N = 41 (3) (1) 
MeV and a^rN ~ 52 (3) (8), respectively. The inclusion of the Nf — 2 a^rN data point greatly reduces the a^^N 
uncertainty as well as brings the two approaches, ABxPT and :^BxPT, closer. In the case of the Nf = 2 + 1 
ensembles, where we fltted solely nucleon mass data, the two approaches give crTrAf-values that differ by 9 MeV. 
This is a novel feature with respect to HBxPT fits where the inclusion of the A-isobar alters the result by more 



For a recent study of the impact of the scale fixing strategy on the value of the sigma terms see Ref. [? 
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than 40 MeV [57|. The inclusion of finite lattice spacing correction to the Nf = 2+1 data tends to reduce 
(T,rjv. Furthermore, we want to call the attention to the fact that our result in Nf = 2 is only cornpatible 
with the experimental determination based on the KA85 nN scattering partial wave analyses of Refs [1^, [SJ] . 
Our Nf = 2 + 1 value is also compatible with the latest determination from the WI08 and EM06 analyses, 
c^ttN = 59(7), which is phenomenologically favored on the grounds of consistency with ttN phenomenology [STj . 
Finally, this Nf = 2 + 1 result would lead, according to the traditional arguments linking sigma terms to the 
baryon-octet mass splittings |7l. Il0l|. to a large strangeness content in the nucleon. However, the uncertainties 
in these arguments have been recently revisited J95| with the conclusion that a a^^N of this size is not at odds 
with, but favored by a negligible strangeness in the nucleon. 

• With both the ABxPT and :;2^BxPT approaches we obtain consistent descriptions of the pion mass dependence 
of the nucleon mass, as can be seen in Figs. [5] and [HI Moreover, for the current IQCD data, all our results are 
compatible within uncertainties and exhibit only small slope variations. However, these small variations translate 
into 11 MeV differences in the value of aj^N at the physical point. For the 2 and 2 + 1 flavor ensembles the Mjr 
distribution of the data points is different. This, together with the differences in the scale-setting procedure 
and the lack of low M^ direct o-TrAr for Nf = 2 + 1 , renders the interpretation of the 1 1 MeV discrepancy based 
on the role of strangeness difficult. To further reduce the uncertainty in the a^^j^ value, IQCD data points with 
smaller uncertainties and less spread would be required. In the Nf = 2 + 1 case a considerable improvement 
could be achieved with a direct measurement of aT^N for A/^ < 300 MeV. It will be interesting to see how the 
Nf = 2 and Nf = 2 + 1 values for (t^^n will change when both data sets become more homogeneous. 
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Appendix A: B^PT Lagrangians 

The counting scheme of Eq. ^ defines the nucleon p'^ self-energy by the sum of the graphs shown in Figs. [T] The 
relevant SU (2) covariant B^PT Lagrangians with explicit A-isobar degrees of freedom are: 
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where the upper indices denote the chiral order. Explicitly, the individual isospin symmetric Lagrangians in 
of external fields and expanded in pion fields tt are: 
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(A8) 
(A9) 



where m^ is the O (p^) pion mass ?7iJ = 2Brri proportional to the chiral condensate B and the current-quark mass 

average m. The Lagrangians C}^'^'' for the nucleon field A^ are those of [9g with a = —4 [8633 + 6115 + eng] a 

combination of C}^1- low energy constants; the £]yJ^ does not produce any nucleon self-energy vertices. The couplings 
of the A-isobar are chosen to be consistent with the covariant construct of the free Rarita-Schwinger theory and hence 
do not contain the unphysical degrees of freedom of vector-spinor fields. The A-isobar Lagrangians and further details 
can be found in |34| - l36l l45l. l97j . There are 13 low energy constants fTro,gAo,ci,C2,C3,C4,cci,HAO:hAO:Mo,AlAO:CiA,ci: 
where C4 and C5 do not contribute to the nucleon mass. 

The loop graphs in Fig. [Ijare divergent in 4 dimensions and need to be regularized. For that we use the dimensional 
regularization with D ~ A — 2e dimensions and renormalize contributions proportional to: 
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For the D-dimensional spin-3/2 propagator we use: 
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The appearing totally anti-symmetric 7-matrices are: 
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Appendix B: Self-energy formulas 

1. Nucleon self-energies 

For the nucleon mass we need the self-energy expressions corresponding to the Feynman-graphs in Fig. [TJ The 
contributions listed in increasing chiral order are: 



S(^) {ml) = SC2 H) , 
S(3)(^2^^) = ^^:,(ml,^)+j:NA3{ml,^) , 
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where we keep the ]i dependence explicit and a 'A' in the index denotes contributions from loop-internal A-isobars. 
The individual unregularized self-energies read: 
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with the expression 
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2. A (1232) self-energies 



In Sec. Ill B I we use the pion mass dependence of the A-isobar to constrain the LEG cia- The A-isobar mass to 
order p^ is 



^A ^ (m^ ) = Mao + Sa2 (ml) + Eaats (™') + SaA3 (ml) , 
where the self-energies are defined as 

K^ {i') = -.9"'' [mi {Mao) + Sf (Mao)] , 
with the unregularized expressions 
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These contributions are the A-isobar versions of the nucleon graphs Ec2, "^NTi and Satas of Fig. [T] The Sca2 is the 
A-isobar contact graph and the Saws and SaA3 are p^ loop with external A-isobars and an internal nucleon and 
A-isobar, respectively. 



3. Finite volume corrections to the nucleon self-energies 

The loop graphs Sats, S7V4, "^ta and Satas, 2]ArA4 of Fig. [T] are subject to finite volume (FV) effects when the 
nucleon is placed in a discretized box. We calculate these effects by the standard techniques of [4J|. In the following 
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we summarize the calculation of the loop-integral with a single propagator and list afterwards all appearing FV 
corrections for the nucleon mass to order p'^. 

For the FV calculation we chose the nucleon rest-frame ^ ~ jopo ~ JoMn- As a consequence all appearing 
loop-integrals can be brought into the form of 
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where no Lorentz-decomposition has to be used, A and B are given 4- vectors and a a power of the 0th- loop momentum 
component. The loop- momentum / is now discretized with respect to the box size L by 
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such that after Wick- rotating and the use of Poisson's formula we get: 
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with j Q Z^ . The case j = corresponds to the usual continuum result whereas the cases j ^ are the finite volume 
corrections. All remaining integrals can be solved analytically. For our nucleon mass expression we need the following 
solutions: 
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where the 1^,^ (x) are modified Bessel- functions of the second kind with F = Ljym? and j = ^/jx+Jy+Jz with 
To collect our final results we use the notations: 
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where the arguments of the self-energies distinguish them from their continuum counterparts. 
The individual finite volume contributions corresponding to the loop-graphs in Fig. [T]are: 
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4. Fit formulas 

In Sees. IIII Al and IIII Bl we use in the x^ fits the foUowmg nucleon mass expressions: 
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+ \aMt + ^C2 {Ml) S^s {Ml) + ^,M^ In ^ 

+EjvA3 {Ml) + ^NAA {Ml) + Ec2 (m') ^'na3 {mI) , (B38) 



where all loops are evaluated at |^ = Mq. The additional terms proportional to ci, as compared to Eq. (J13p . come 
from the discussion in Sec. IIIBI In the case of fits with finite volume corrections, we add the following expressions: 
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